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ABSTRACT 


An  experimental  investigation  of  heat  transfer  to  hemispheres  with 
both  catalytic  and  noncatalytic  surfaces  in  nonequilibrium  dissociated 
hypersonic  nitrogen  flow  has  been  conducted.  The  wind  tunnel  flow  con¬ 
ditions  provided  low  Reynolds  numbers,  allowing  the  assumption  of 
completely  frozen  shock  and  boundary  layers  while  introducing  the  in¬ 
fluences  of  second-order  boundary  layer  effects  such  as  vorticity  inter¬ 
action.  The  results  for  copper  surfaces,  when  correlated  using  the 
stagnation  point  heat  transfer  equation  of  Rosner  together  with  Lees' 
heat  transfer  distribution,  are  in  good  agreement  with  the  second-order 
theory  of  Cheng  under  the  assumption  of  complete  surface  recombination 
of  atoms.  Reductions  in  heat  transfer  rate  up  to  30  percent  were  obtained 
with  the  use  of  noncatalytic  coatings  applied  to  the  probes. 
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SECTION  I 
INTRODUCTION 


It  has  been  irecognized  for  some  time  that  a  body  traveling  at  hyper¬ 
sonic  speeds  through  the  upper  atmosphere  may  experience  a  mechanism 
of  heat  transfer  other  than  the  purely  molecular  heat  conduction  process 
of  conventional  fluid  mechanics.  This  will  occur  when  the  shock  wave 
ahead  of  a  body  converts  a  portion  of  the  flight  kinetic  energy  into  the 
chemical  energy  of  dissociation  of  the  air  molecules.  In  this  event,  two 
driving  mechanisms  for  heat  transfer  to  the  body  become  apparent;  the 
temperature  gradient  across  the  boundary  layer  and  the  concentration 
gradient  of  atoms  within  the  boundary  layer.  The  relative  importance  of 
these  two  mechanisms  in  fixing  the  amount  of  heat  transfer  which  will 
occur  is  determined  by  the  conditions  within  the  boundary  layer  and  at 
the  body  surface. 

If  the  boundary  layer  flow  properties,  particularly  the  density,  are 
such  that  the  characteristic  time  required  for  atom  recombination  is  very 
much  smaller  than  the  time  required  for  atom  diffusion  across  the  bound¬ 
ary  layer,  then  the  equilibrium  boundary  layer  exists,  in  which  the  re¬ 
combination  is  completed  before  the  atoms  can  diffuse  to  the  cold  surface. 
In  the  other  extreme,  the  "frozen"  boundary  layer  case,  the  characteristic 
time  for  atom  recombination  is  so  large  that  no  recombination  can  occur 
before  the  atoms  have  diffused  to  the  surface.  In  these  regimes,  and  in 
the  intermediate  regime  where  the  characteristic  times  for  diffusion  and 
recombination  are  of  comparable  order  of  magnitude,  the  dimensionless 
combination  of  properties  known  as  the  Lewis  number  is  of  particular 
importance. 

The  dimensionless  group  p  CpD/k,  commonly  called  the  Lewis  num¬ 
ber,  arises  as  a  matter  of  convenience  in  the  analysis  of  convective  heat 
transfer  from  dissociated  gases.  It  represents  the  ratio  of  the  diffusivity 
of  dissociated  atoms  to  the  thermal  diffusivity  of  heat  by  conduction.  For 
similar  profiles  of  temperature  and  atom  concentration  across  the  bound¬ 
ary  layer,  and  for  a  given  total  enthalpy  at  the  edge  of  the  boundary  layer, 
the  heat  transfer  rate  is  strictly  dependent  upon  the  enthalpy  difference 
when  Le  =  1,  regardless  of  the  proportion  of  the  enthalpy  gradient  asso¬ 
ciated  with  a  temperature  gradient  or  an  atom  concentration  gradient. 

For,  with  Le  =  1,  energy  will  be  transported  across  the  boundary  layer  by 
both  conduction  and  atom  diffusion  at  the  same  rate  per  unit  temperature 
gradient  or  concentration  gradient.  It  will  be  observed  that  in  the  restric¬ 
tion  of  similar  profiles  of  temperature  and  atom  concentration,  the  condi¬ 
tion  of  zero  atom  concentration  at  the  solid  surface  is  analagous  to  the  zero 
temperature  jump  condition. 
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For  the  more  general  case,  when  the  Lewis  number  is  different  from 
unity,  the  heat  transfer  rate  is  no  longer  simply  dependent  upon  the  total 
enthalpy  difference  across  the  boundary  layer.  The  greater  atom  dif¬ 
fusion  rate  associated  with  a  Lewis  number  greater  than  unity  augments 
the  transport  of  energy  to  the  cold  surface,  whereas  the  smaller  atom 
diffusion  rate  associated  with  a  Lewis  number  less  than  unity  retards  the 
transport  of  energy  to  the  surface.  For  most  partially  dissociated 
gases  of  aerodynamic  interest,  including  air,  the  Lewis  number  is 
slightly  greater  than  unity,  the  effect  being  to  increase  the  heat  trans¬ 
fer  rate,  say,  on  the  order  of  10  percent. 

The  preceding  statements  have  been  confined  strictly  to  the  case  in¬ 
volving  similar  temperature  and  atom  concentration  profiles;  that  is, 
the  atom  concentration  is  zero  at  the  cold  surface.  For  a  partially  or 
completely  frozen  boundary  layer,  the  case  of  zero  atom  concentration 
at  the  wall  represents  the  limiting  case  of  a  surface  which  is  fully  cata¬ 
lytic  to  atom  recombination;  that  is,  all  atoms  which  diffuse  to  the  sur¬ 
face  recombine  there,  depositing  their  chemical  energy  on  the  surface. 

In  general,  a  surface  will  not  be  fully  catalytic.  Therefore,  some  of  the 
atoms  which  reach  the  surface  will  not  recombine  there,  resulting  in  a 
finite  atom  concentration  at  the  surface.  The  obvious  result,  if  the 
boundary  layer  is  partially  or  completely  frozen,  is  that  a  surface  which 
is  not  fully  catalytic  will  cause  the  heat  transfer  rate  to  be  reduced  be¬ 
cause  of  the  decrease  in  surface  recombination  and  the  resulting  decrease 
in  the  atom  concentration  gradient.  The  condition  of  completely  frozen 
boundary  layer  and  completely  non- catalytic  surface  represents  the  lower 
limit  on  the  heat  transfer,  where  the  chemical  contribution  is  zero. 

The  laminar  heat  transfer  to  the  stagnation  region  of  a  blunt-nosed 
body  in  hypersonic  flow  has  been  treated  theoretically  by  several  in¬ 
vestigators  (Refs.  1  through  6).  Lees  (Ref.  1)  considered  only  the  equi¬ 
librium  boundary  layer  and  the  completely  frozen  boundary  layer  with  a 
fully  catalytic  surface.  Approximate,  closed-form  solutions  for  these 
two  cases  were  obtained  from  the  boundary  layer  equations  simplified  on 
the  basis  of  physical  arguments.  Fay  and  Riddell  (Ref.  2)  obtained  nu¬ 
merical  solutions  to  the  boundary  layer  equations  over  homogeneous  (gas- 
phase)  recombination  rates  from  equilibrium  flow  to  frozen  flow  in  the 
boundary  layer.  Again,  these  results  apply  only  for  the  limiting  case  of 
a  fully  catalytic  surface. 

Goulard  (Ref.  3)  integrated  the  frozen  laminar  boundary  layer  equa¬ 
tions  and  obtained  a  solution  to  the  stagnation  heat  transfer  problem  with 
an  arbitrary  degree  of  catalytic  activity  at  the  surface.  This  was  done 
simply  by  introducing  a  correction  factor,  <f>,  which  is  a  function  of  the 
flow  conditions,  nose  geometry,  and  wall  catalytic  reaction  rate  constant. 
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to  account  for  the  realistic  condition  in  which  the  atom  concentration  at 
the  wall  is  not  zero.  Rosner  (Ref.  4)  has  suggested  an  expression  for 
the  frozen  boundary  layer  stagnation  heat  transfer  which  is  quite  simi¬ 
lar  to  the  results  of  both  Goulard  and  Fay  and  Riddell,  while  employing 
a  more  consistent  use  of  the  correction  factor,  6,  when  the  surface 
temperature  cannot  be  neglected.  This  expression  is  considered  in  de¬ 
tail  in  a  later  section  of  this  report. 

Several  other  authors  have  treated  analytically  the  problem  of 
stagnation  point  heat  transfer  in  hypersonic  flow.  Chung  and  Liu  (Ref.  5) 
have  developed  an  approximate  analysis  to  predict  the  heat  transfer  to  a 
noncatalytic  surface.  Their  results  were  then  generalized  to  apply  to 
simultaneous  gas-phase  recombination  and  surface  catalytic  recombina¬ 
tion.  Inger  (Ref.  8)  has  also  presented  an  approximate  theory  of  non¬ 
equilibrium  stagnation  point  boundary  layers  with  atom  recombination 
at  the  surface. 

Numerous  other  authors  (Ref.  7)  have  investigated  various  facets  of 
the  effects  of  chemical  kinetics  in  convective  heat  transfer.  Nearly  all 
of  these  investigations  have  been  of  a  theoretical  nature,  and  the  prin¬ 
ciples  involved  in  the  theory  are  well  understood.  In  fact,  Spaulding 
(Ref.  8)  has  stated  that  "no  mysteries  now  remain"  in  our  ability  to  pre¬ 
dict  the  effects  of  chemical  reaction  on  heat  transfer  from  gases. 

On  the  other  hand,  the  view  taken  by  the  present  writer  is  that  there 
has  been  far  too  little  experimental  verification  of  the  existing  theories 
on  heat  transfer  from  chemically  reacting  gases.  This  lack  of  verifica¬ 
tion  is  not  surprising  in  view  of  the  difficulties  involved  in  making  such 
measurements. 

Rose  and  Stark  (Ref.  9)  have  reported  the  results  of  shock  tube  meas¬ 
urements  of  stagnation  point  heat  transfer  in  dissociated  air.  They 
obtained  good  agreement  with  the  equilibrium  boundary  layer  theory  of 
Fay  and  Riddell  (Ref.  2).  The  flow  conditions  were  such  that  any  non¬ 
equilibrium  effects  in  the  boundary  layer  and  at  the  surface,  if  present 
at  all,  were  very  small.  Busing  (Ref.  10)  obtained  shock  tube  measure¬ 
ments  of  the  stagnation  point  heat  transfer  to  both  catalytic  and  noncata¬ 
lytic  walls.  Although  his  results  showed  a  reduction  in  heat  transfer  with 
the  noncatalytic  surface,  the  size  of  the  reduction  was  not  nearly  as  large 
as  expected. 

The  differential  heat  transfer  to  catalytic  and  noncatalytic  surfaces 
in  a  steady- state,  supersonic  wind  tunnel  was  measured  by  Rosner 
(Ref.  11)  and  by  Winkler  and  Griffin  (Ref.  12).  A  sizeable  heat  transfer 
difference  was  obtained  in  each  investigation.  These  experiments, 
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however,  were  not  designed  to  allow  comparison  of  the  results  with  an 
appropriate  theory  for  heat  transfer. 

Hartunian  and  Thompson  (Ref.  13)  have  reported  preliminary  meas¬ 
urements  of  catalytic  and  noncatalytic  differential  heat  transfer  in  a  low 
density  shock  tube.  Their  large  (17-in.  -diam)  shock  tube  permits  larger 
probes  and  greater  test  times  than  are  usually  possible  with  smaller, 
low  pressure  shock  tubes. 

The  need  for  additional  experimental  data  on  nonequilibrium  heat 
transfer  with  surface  catalysis  is  clear.  The  available  experimental 
results  generally  have  not  been  compared  with  an  appropriate  theory 
(such  as  Fay-Riddell,  Ref.  2)  for  predicting  the  absolute  level  of  heat  ■ 
transfer.  The  usual  procedure  has  been  simply  to  show  the  percentage 
reduction  in  heat  transfer  caused  by  a  noncatalytic  wall,  with  little  con¬ 
cern  for  the  absolute  values.  It  is  fair  to  say  that  most  experimental 
results  of  nonequilibrium  heat  transfer  to  catalytic  and  noncatalytic  sur¬ 
faces  have  been  qualitative,  not  quantitative. 

The  purpose  of  this  report  is  to  present  quantitative  results  of  an 
experimental  investigation  of  the  heat  transfer  to  blunt  bodies,  with  both 
catalytic  tnd  noncatalytic  surfaces,  in  nonequilibrium  dissociated  hyper¬ 
sonic  flow.  These  experiments  were  performed  in  a  small,  continuous 
flow,  low  density  hypersonic  wind  tunnel  in  operation  at  the  von  Karman 
Gas  Dynamics  Facility  (VKF). 


SECTION  II 
FLOW  CONDITIONS 


2.1  THE  WIND  TUNNEL 

The  wind  tunnel  used  in  this  investigation  consists  of  the  following 
major  components,  in  stream  wise  order: 

1.  Direct  current  arc -heater  (Thermal  Dynamics  U-50),  with  40-kw 
power  supply; 

2.  Settling  section  of  variable  size  but  normally  of  3 -in.  diam  and 
6-  to  10-in.  length; 

3.  Aerodynamic  nozzle  of  variable  size  with  0.  10-  to  0.  75-in.  -diam 
throat  and  2.0-  to  8.  0-in.  -diam  exit; 

4.  A  tank  of  48- in.  diam  surrounding  the  test  section  and  containing 
instrumentation  and  probe  carrier; 

5.  Interchangeable  diffuser; 

6.  Water-cooled  heat  exchanger; 
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7.  Air- ejector  of  two  stages;  and 

8.  The  VKF  mechanical  vacuum-pumping  system. 

A  photograph  of  the-tunnel  is  presented  as  Fig.  1. 

Since  this  wind  tunnel  is  used  to  study  aerodynamic  problems  of  low 
density,  hypersonic  flows,  it  is  necessary  that  the  tunnel  flow  conditions 
be  known  with  a  higher  degree  of  accuracy  than  is  common  with  arc- 
heated  facilities  normally  used  in  heat  transfer  experiments.  A  survey 
of  some  of  the  techniques  used  in  flow  calibration  and  diagnosis  is  pre¬ 
sented  in  Refs.  14  and  15.  These  techniques  have  been  used  in  the 
present  investigation  to  establish  the  test  section  flow  conditions  with 
what  is  believed  to  be  a  high  degree  of  accuracy. 

The  wind  tunnel  is  usually  operated  with  nitrogen  at  reservoir  pres¬ 
sures  greater  than  one  atmosphere  and  at  a  test  section  total  enthalpy 
corresponding  to  an  equilibrium  temperature  of  about  3000°K.  It  has 
been  shown  (Ref.  15)  that  the  gas,  at  this  pressure  and  enthalpy,  will 
approach  thermodynamic  equilibrium  upstream  of  the  nozzle  throat  pro¬ 
vided  the  settling  chamber  is  of  sufficient  length  to  allow  complete 
mixing  and  recombination  of  the  highly  energetic  species  produced  by 
the  arc  heater,  where  temperatures  in  excess  of  those  required  for  nitro¬ 
gen  dissociation  are  present.  Therefore,  the  number  of  nitrogen  atoms 
in  the  test  section  of  the  nozzle  is  negligible  at  the  usual  tunnel  operating 
conditions.  However,  when  the  tunnel  is  operated  at  a  reservoir  pres¬ 
sure  much  less  than  one  atmosphere,  there  is  an  insufficient  number  of 
collisions  within  the  reservoir  to  allow  complete  atom  recombination 
upstream  of  the  nozzle  throat.  Hence,  a  small  fraction  of  the  test  section 
flow  will  be  dissociated  even  when  the  average  temperature  is  well  below 
the  dissociation  level,  resulting  in  a  significant  increase  in  the  total 
enthalpy  above  the  equilibrium  value  (cf,  Ref.  15).  It  is  this  capability, 
together  with  the  refined  calibration  procedures  and  the  features  of  con¬ 
tinuous  operation,  which  makes  this  facility  attractive  for  the  study  of 
heat  transfer  in  dissociated  nonequilibrium  flow. 

2.2  FLOW  DIAGNOSIS 

A  wide  variety  of  flow  conditions  involving  nonequilibrium  dis¬ 
sociation  may  be  obtained  in  Tunnel  L  by  proper  choice  of  the  nitrogen 
flow  rate,  total  pressure,  settling  chamber  geometry,  and  the  nozzle 
throat  area.  The  feasibility  of  the  present  investigation  was  demon¬ 
strated  in  a  brief  preliminary  investigation*  by  the  present  author  and 


♦These  preliminary  results  are  reported  as  an  example  in  Ref.  16. 
They  are  qualitatively  correct;  however,  the  numerical  values  may  be 
slightly  revised  in  light  of  the  present  experiments. 
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a  colleague,  J.  T.  Miller.  The  same  flow  conditions  used  in  this  early 
investigation  were  selected  for  the  present  investigation.  The  nitrogen 
flow  rate  was  15.  2  lbm/hr  at  a  measured  reservoir  pressure  of 
3.  5  psia.  A  settling  chamber  with  a  diameter  of  3  in.  and  nominal 
length  of  8  in.  was  used.  The  aerodynamic  nozzle  had  a  throat  diameter 
of  0.  397  in.  and  a  conical  expansion  section  with  a  15-deg  half-angle. 

2.2.1  Total  Enthalpy  Measurements 

When  the  thermo-chemical  flow  conditions  in  a  wind  tunnel  nozzle 
are  not  known  a  priori,  it  is  then  impossible  to  compute  the  total  en¬ 
thalpy  of  the  flow  from  the  usual  measured  quantities  such  as  mass  flow 
rate,  total  pressure,  and  throat  area.  The  investigator  must  turn  to  a 
direct  measurement  of  the  total  enthalpy  in  the  test  section. 

The  device  used  for  this  measurement  in  the  present  investigation 
is  the  total -enthalpy  and  mass-flow  probe  shown  in  Fig.  2.  The  probe 
consists  of  two  concentric  water-cooled  jackets,  separated  by  an  air 
space  to  provide  thermal  insulation.  Each  jacket  is  knife- edged  at  the 
front  to  minimize  the  contact  area  for  heat  conduction  between  jackets 
and  to  provide  the  proper  geometry  to  enable  the  bow  shock  to  be 
"swallowed"  when  the  probe  is  placed  axially  in  the  test  section.  The 
function  of  the  inside  jacket  is  to  remove  energy  from  the  captured 
flow  and  deposit  this  energy  in  the  cooling  water,  whose  flow  rate  and 
temperature  rise  are  measured.  The  total  energy  convected  into  the 
calorimeter  is  then  the  sum  of  the  energy  increase  in  the  cooling  water 
of  the  internal  jacket  and  the  residual  energy  of  the  gas  leaving  the 
calorimeter,  as  measured  by  a  thermocouple. 

The  total  enthalpy  of  the  nozzle  flow  is  obtained  by  dividing  the  total 
energy  convected  into  the  calorimeter  by  the  mass  flow  captured  by  the 
calorimeter.  If  the  bow  shock  is  completely  swallowed,  the  mass  flow 
rate  is  simply  the  quantity  times  the  probe  capture  area.  If  the 

bow  shock  is  not  swallowed,  the  captured’ mass  flow  rate  must  be  meas¬ 
ured  independently. 

This  probe  was  designed  to  allow  an  accurate  measurement  of  the 
captured  flow  rate.  The  flow  is  directed  from  the  probe  to  an  evacuated 
vessel  of  known  volume.  The  flow  rate  is  determined  by  observing  the 
pressure  rise  in  the  vessel  during  a  measured  time  interval.  The  de¬ 
tails  of  this  procedure  have  been  reported  in  Ref.  14.  The  measured 
flow  rate  was  found  to  be  unchanged  when  the  time  interval  and  the  back 
pressure  were  varied  over  a  wide  range,  indicating  that  the  bow  shock  at 
the  probe  entrance  was  indeed  swallowed. 
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Since  it  is  virtually  impossible  to  eliminate  all  internal  paths  for 
heat  loss  in  a  probe  of  this  size  in  such  an  extreme  environment,  this 
device  is.  an  imperfect  calorimeter  and  must  be  calibrated.  The  cali¬ 
bration  procedure  consisted  of  testing  the  probe  in  the  flows  from  two 
other  nozzles,  operating  at  three  flow  conditions  for  which  the  total 
enthalpy  is  well  known  from  other  measurements.  The  probe  measure¬ 
ments  were  found  to  be  7.  8,  5,  3,  and  9.  1  percent  higher  than  the  actual 
total  enthalpy  for  the  three  flow  conditions.  The  calibration  correction 
factor  used  for  the  nonequilibrium  flow  condition  was  7.  4  percent,  the. 
average  of  the  three  deviations.  Using  this  factor  to  adjust  the  probe 
measurement,  the  total  enthalpy  of  the  dissociated  flow  was  found  to  be 
1130  Btu/lbm. 

If  it  is  further  assumed  that  all  of  the  total  enthalpy  tied  up  in  non¬ 
equilibrium  modes  at  the  nozzle  throat  is  caused  by  dissociation,  i.  e. , 
vibration  is  in  equilibrium  at  the  nozzle  throat,  then  the  chemical 
energy  of  dissociation  is  computed  to  be  393  Btu/lbm.  Since  the  dis¬ 
sociation  energy  for  nitrogen  is  14,  270  Btu/lbm.  the  dissociation  frac¬ 
tion,  «,  is  found  to  be  0.  0275. 

2.2.2  Nitric  Oxide  Titration 

In  order  to  obtain  an  independent  measurement  of  the  nitrogen  atom 
fraction  in  the  test  section,  the  nitric  oxide  light  titration  technique  was 
used.  This  technique,  first  described  by  Spealman  and  Rodebush 
(Ref.  17)  and  discussed  later  in  more  detail  by  Kistiakowsky  and  Volpi 
(Ref.  18),  is  based  on  the  fact  that  the  reaction 

N  +  NO  -  N2  +  0 

is  very  fast  compared  with  the  reaction  of  nitrogen  atoms  with  any  other 
particles.  The  addition  of  nitric  oxide  to  a  stream  containing  nitrogen 
atoms  is  accompanied  by  a  pronounced  change  in  the  emission  spectrum, 
producing  a  change  of  colors  from  orange-yellow  through  pink  and 
purplish-blue  to  colorless  and  then  to  whitish-green  as  increasing  amounts 
of  nitric  oxide  are  added  (Refs.  19  through  22).  The  titration  end  point, 
or  the  condition  at  which  all  of  the  nitric  oxide  is  consumed  by  the  nitrogen 
atoms,  is  just  reached  when  the  emission  of  light  disappears.  The  titra¬ 
tion  end  point  is  very  sharply  defined,  since  only  a  small  increase  in  the 
nitric  oxide  flow  rate  will  result  in  the  whitish- green  emission  of  the 
reaction 

NO  +  O-NO,  +  hv 

which  occurs  whenever  nitric  oxide  is  added  in  excess  of  the  stoichiomet¬ 
ric  amount  for  reaction  with  nitrogen  atoms.  The  mole  flow  rate  of 
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nitrogen  atoms  is  determined  by  measuring  the  nitric  oxide  flow  rate 
at  the  titration  end  point,  since  the  reaction  occurs  in  equimolar  pro¬ 
portions.  Discussions  of  the  merits  of  the  nitric  oxide  titration 
technique  may  be  found  in  Refs.  23  and  24. 

The  nitric  oxide  was  injected  into  the  test  section  through  a  stain¬ 
less  steel,  uncooled,  1/4-in.  -diam  probe  positioned  on  the  axis  of 
symmetry.  The  upstream  end  of  the  probe  was  closed.  Several  holes 
were  drilled  in  the  probe  at  a  position  about  1/2  in.  behind  the  tip  to 
permit  a  uniform  radial  injection  of  the  nitric  oxide. 

It  was  found  that  the  emission  of  light  from  the  tank  volume  sur¬ 
rounding  the  high  speed  portion  of  the  flow  was  very  sensitive  to  the 
nitric  oxide  flow  rate,  while  the  colors  produced  in  the  high  speed  flow 
of  nitrogen  varied  from  point  to  point  in  the  stream.  This  color  varia¬ 
tion  undoubtedly  resulted  from  the  poor  mixing  obtained  in  the  immediate 
vicinity  of  the  injection  probe,  complicated  by  the  possibility  that  the 
chemical  reactions  which  produce  the  afterglow,  even  though  very  fast, 
were  not  completed  near  the  probe.  The  emission  of  light  in  all  parts  of 
the  tank  volume  outside  the  high  speed  stream  indicated  that  a  significant 
quantity  of  nitrogen  atoms  and  nitric  oxide  molecules  was  reaching  this 
area.  It  has  been  shown  (Ref.  25)  that  there  is  a  reverse  flow  along  the 
walls  of  the  convergent  diffuser  entrance  section.  This  reverse  flow  has 
been  observed  by  attaching  tufts  to  the  diffuser  walls.  It  was  also  plainly 
visible  during  the  whitish- green  emission  caused  by  excessive  nitric 
oxide.  Therefore,  this  reverse  flow  and  mixing  phenomenon  is  the 
mechanism  by  which  the  reacting  species  reach  the  outer  volume  of  the 
tank. 

With  no  nitric  oxide  injection  into  the  stream,  the  yellow-orange 
Lewis -Rayleigh  nitrogen  afterglow  was  clearly  visible  when  the  room 
lights  were  darkened.  This  emission,  which  accompanies  the  homo¬ 
geneous  decay  of  nitrogen  atoms  (Refs.  22  and  26),  filled  the  entire  tank 
surrounding  the  test  section.  As  a  small  amount  of  nitric  oxide  was  in¬ 
jected  into  the  flow,  the  emission  inside  the  tank  became  darker  in  color 
and  less  intense  until,  at  a  particular  nitric  oxide  flow  rate,  the  visible 
emission  was  completely  extinguished.  A  very  minute  increase  in  the 
nitric  oxide  flow  rate  resulted  in  a  series  of  whitish- green  flashes  of 
light  which  traversed  the  entire  volume  of  the  tank  in  a  random  manner. 

A  further  increase  in  the  nitric  oxide  flow  rate  caused  the  emission  of 
the  relatively  intense  whitish- green  light  continuously  throughout  the 
entire  tank  volume.  The  nitric  oxide  flow  rate  was  measured  at  the  con¬ 
dition  where  the  light  was  extinguished,  just  prior  to  the  appearance  of 
the  white  flashes. 
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The  results  of  the  titration  procedure,  in  terms  of  the  measured 
nitrogen  atom  fraction,  are  presented  in  Fig.  3  for  various  probe  posi¬ 
tions  in  the  test  section.  The  upstream  tip  of  the  probe,  located  on  the 
flow  centerline,  was  moved  axially  to  several  stations  in  the  test  sec¬ 
tion  near  the  nozzle  exit.  The  measured  nitrogen  atom  fraction  decreased 
as  the  probe  was  moved  downstream  from  the  nozzle  exit.  Although  this 
result  is  in  the  correct  direction  to  account  for  atom  recombination  in  the 
nozzle  flow,  the  better  explanation  is  found  in  the  mixing  process  down¬ 
stream  of  the  probe  and  at  the  diffuser.  Since  the  titration  end  point  is 
determined  by  observing  the  color  changes  in  the  tank  volume  outside  the 
primary  flow  region,  the  technique  will  be  accurate  only  if  the  nitric  oxide 
and  nitrogen  atoms  and  their  products  of  reaction  are  well  mixed  at  the 
diffuser  entrance  and  are  trapped  in  the  reverse  flow  at  the  diffuser  walls 
in  the, proper  relative  proportion.  By  injecting  the  nitric  oxide  at  the  tank 
wall  instead  of  on  the  flow  centerline,  it  was  found  that  far  less  titrant  gas 
was  required  to  reach  the  apparent  titration  point.  Therefore,  since  the 
overwhelming  percentages  of  nitrogen  atoms  and  titrant  gas  pass  through 
the  diffuser  without  being  recirculated,  the  results  of  the  technique  are  de¬ 
pendent  upon  the  efficiency  of  the  mixing  process.  The  change  in  the  probe 
location  undoubtedly  affects  the  mixing  process  and  hence  the  measured 
atom  fractions.  Even  though  the  titration  technique  must  be  considered 
only  semi-quantitative  in  the  present  application,  its  results  were  very  use¬ 
ful  in  increasing  confidence  in  the  atom  fraction  measurement  obtained  with 
the  calorimeter  probe.  The  calorimeter  probe  measurement  of  0.  0275 
falls  amazingly  close  to  the  center  of  the  band  of  measurements  shown  in 
Fig.  3. 

2.2.3  Impact  Pressure  Survey 

The  impact  pressure  probe  is  of  major  importance  in  flow  calibration 
since  the  flow  conditions  can,  in  principle,  be  calculated  from  the  meas¬ 
ured  impact  and  total  pressures  if  the  correct  thermo-chemical  flow  model 
is  known. 

In  this  investigation  the  test  section  flow  properties  (density,  velocity, 
and  Mach  number)  were  calculated  from  the  results  of  an  impact  pressure 
survey  in  the  test  section.  For  this  calculation  it  was  assumed  that  the 
previously  measured  value  of  the  dissociation  fraction  («  =  0.  0275)  re¬ 
mained  constant  throughout  the  entire  nozzle,  while  the  molecular  vibration 
mode  was  in  equilibrium  upstream  of  the  throat  and  frozen  downstream  of 
the  throat.  This  assumption  is  at  least  qualitatively  correct  since  the  num¬ 
ber  of  collisions  required  for  atom  recombination  is  several  orders  of 
magnitude  greater  than  the  number  required  for  vibrational  relaxation. 

Asa  further  check  on  the  validity  of  the  assumed  thermo- chemical 
flow  model,  use  was  made  of  a  solution  recently  developed  by  Kinslow  and 
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Miller,  colleagues  of  the  writer.  This  solution  (Ref.  27)  applies  the 
basic  equations  for  the  conservation  of  mass,  momentum,  and  energy 
to  the  one- dimensional  adiabatic  flow  of  an  inviscid  diatomic  gas  with 
finite  vibrational  relaxation  and  dissociation  recombination  rates. 
Together  with  the  rate  equation,  these  equations  are  solved  by  an 
iteration  technique  on  a  digital  computer  from  any  initial  condition  for 
the  flow  through  a  hypersonic  nozzle.  For  the  conditions  measured  in 
the  test  section  in  this  investigation  (h0  =  1130  Btu  lbm,  a  =  0.  0275), 
it  was  determined  that  the  dissociation  was  indeed  frozen  throughout 
the  nozzle,  and  that  molecular  vibration,  while  frozen  downstream  of 
the  nozzle  throat,  was  not  too  different  from  its  equilibrium  value 
upstream  of  the  throat.  Therefore,  the  assumed  flow  model  was  verified 
to  be  substantially  correct,  validating  the  flow  property  calculation  pro¬ 
cedure  based  on  the  measured  impact  pressures. 

Some  of  the  important  flow  parameters  are  given  in  Table  I  for 
several  locations  along  the  flow  centerline. 

2.3  BOUNDARY  LAYER  CONDITIONS 
2.3.1  Ra-Entry  Simulation 

The  dissociation  which  may  exist  within  the  boundary  layer  around 
an  actual  re-entry  vehicle  is  generated  by  the  detached  shock  wave  as 
it  passes  through  the  ambient  atmosphere.  If  the  flow  field  behind  the 
shock  wave  is  in  equilibrium,  then  the  highest  possible  amount  of  dis¬ 
sociation  will  have  been  produced.  If  the  boundary  layer  flow  remains 
in  equilibrium,  then  the  computation  of  the  heat  transfer  for  this  case 
is  straightforward  since,  for  the  cold-wall  condition  considered  here, 
the  effects  of  surface  atom  recombination  are  absent. 

The  more  interesting  (and  difficult)  situation  occurs  when  the  flow 
field  behind  the  shock  wave  is  not  in  equilibrium.  In  this  case,  the  shock- 
induced  dissociation  will  lag  behind  the  equilibrium  value  up  to  the  outer 
edge  of  the  boundary  layer.  Then,  if  the  boundary  layer  itself  is  not  in 
equilibrium,  at  some  point  within  the  boundary  layer  the  actual  dissocia¬ 
tion  fraction  will  reach  a  maximum  value  and  then  decrease  toward  the 
wall.  The  computation  of  the  heat  transfer  for  this  case  is  difficult 
because,  in  addition  to  the  effects  of  wall  recombination,  the  quantity 
of  energy  deposited  in  dissociation  is  unknown. 

Exact  duplication  of  the  actual  thermo- chemical  re-entry  model  in 
a  wind  tunnel  is  difficult  to  achieve.  In  addition  to  the  difficulties  associ¬ 
ated  with  obtaining  the  high  temperatures  necessary  for  shock-induced 
dissociation,  the  low  densities  required  for  a  nonequilibrium  boundary 
layer  would  result  in  a  nonequilibrium  nozzle  expansion  and  a  significant 
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amount  of  free-stream  dissociation.  Furthermore,  the  difficulties  in 
experimentally  determining  the  actual  thermo- chemical  state  of  the  gas 
near  the  body  appear  formidable. 

The  alternative  procedure  for  the  wind  tunnel  simulation  of  re-entry 
is  the  production  of  a  nonequilibrium  dissociated  free-stream  flow,  as 
was  done  in  the  present  work.  The  total  temperature  corresponding  to 
the  translational,  rotational,  and  vibrational  degrees  of  freedom  was 
about  1525  °K.  Obviously  there  could  be  no  additional  dissociation  pro¬ 
duced  through  the  shock  region  at  this  temperature,  but  the  possibility 
of  recombination  in  this  region  must  be  examined. 

The  question  of  gas -phase  atom  recombination  in  the  stagnation 
region  of  blunt  bodies  at  low  Reynolds  numbers  has  been  considered  by 
Goulard  (Ref.  3),  Grier  and  Sands  (Ref.  28),  Rosner  (Ref.  29),  Cheng  (30), 
and  Buckmaster  (Ref.  31).  Cheng  assumed  a  single  dissociation- 
recombination  reaction  of  the  type  N2  +  M  t  N+N+M,  where  M  represents 
a  third  body.  He  then  calculated  the  value  of  the  dissociation  fraction 
through  the  shock  layer  for  several  values  of  Reynolds  number.  At  the 
lowest  Reynolds  number  considered,  he  found  that  both  dissociation  and 
recombination  were  negligible  in  the  shock  layer.  As  the  degree  of  rare¬ 
faction  decreased,  the  effect  of  dissociation  in  the  outer  portion  of  the 
shock  layer  became  significant  while  recombination  remained  negligible. 
As  pointed  out  by  Buckmaster,  this  is  primarily  because  recombination 
requires  three-body  collisions,  which  occur  less  frequently  than  the 
dissociation  collision.  At  the  highest  Reynolds  number  considered  by 
Cheng,  which  was  higher  than  those  of  the  present  investigation,  the 
gas-phase  recombination  was  still  negligible.  Comparison  of  the  present 
flow  conditions  with  the  numerical  results  of  Rosner  (Ref.  29)  for  a 
similar  flow  condition  results  in  the  same  conclusion,  that  is,  negligible 
recombination  in  the  boundary  layer.  Furthermore,  as  Rosner  points 
out,  if  recombination  is  not  appreciable  in  the  boundary  layer  itself, 
it  will  certainly  be  negligible  in  the  entire  shock  layer  region, 

2.3.2  S»cond>Ord*r  Effects 

In  the  low  density  regime  in  which  the  present  experiments  were 
conducted,  the  usual  first-order  boundary  layer  theory  must  be  modified 
to  consider  second-order  effects  such  as  external  vorticity,  displace¬ 
ment,  curvature,  slip,  and  temperature  jump.  A  number  of  investigators 
■(Refs.  32  through  45)  have  considered  one  or  more  of  these  second-order 
effects  using  various  techniques.  Van  Dyke  (Ref.  46)  has  presented  an 
excellent  review  and  critical  discussion  of  this  subject.  It  is  sufficient 
to  state  that  the  agreement  between  the  several  theories  has  not  always 
been  good.  Similarly,  the  agreement  between  experiment  and  theory  is 
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poor  in  many  instances.  Furthermore,  there  is  some  disagreement 
between  experiments  performed  at  different  laboratories. 

It  is  hoped  ihat  the  experiments  reported  here  will  be  welcomed 
not  only  by  those  whose  main  interests  lie  in  the  chemical-kinetic  and 
surface  effects  of  heat  transfer  but  also  by  those  who  are  concerned 
with  the  second-order  effects  in  hypersonic  boundary  layers.  Indeed, 
in  many  physical  situations  these  two  classes  of  problems  will  appear 
simultaneously. 


SECTION  III 

HEAT  TRANSFER  PROBES 

The  total  heat  transfer  rate  to  the  hemispherical  noses  of  water- 
cooled,  calorimeter -type  probes  was  measured.  This  type  probe, 
shown  in  Fig.  4,  measures  the  steady-state,  linear  temperature  drop 
along  a  conductor  of  known  thermal  conductivity  and  cross-sectional 
area.  The  steady-state  heat  transfer  rate  to  the  hemispherical  nose  is 
then  calculated  using  Fourier's  law.  Since  the  conductor  within  the  probe 
is  shielded  from  the  external  flow,  the  heat  losses  are  negligible.  The 
probe  is  calibrated  by  applying  a  known  heat  input,  and  a  calibration 
factor  of  unity  is  consistently  obtained  at  the  author's  laboratory  for 
this  type  probe. 

In  order  to  obtain  data  over  a  wider  range  of  Reynolds  numbers, 
three  probes,  with  nose  radii  of  0.  25,  0.50,  and  0.75  in.,  were  used. 
The  hemispherical  nose  of  each  probe  was  machined  from  solid  copper 
of  high  purity.  This  copper  surface  was  used  as  the  catalytic  probe, 
with  no  surface  preparation  other  than  polishing  with  fine  emery  paper. 

The  same  three  probes  were  used  to  obtain  the  noncatalytic  wall 
data.  Two  noncatalytic  surfaces  were  tested.  The  first  noncatalytic 
surface  was  ortho-phosphoric  acid  (H3PO4),  applied  to  the  copper  sur¬ 
face  with  a  brush,  and  allowed  to  dry  under  vacuum  conditions.  The 
second  noncatalytic  surface  was  a  layer  of  silicon  monoxide,  evaporated 
onto  the  copper  hemisphere.  a  The  thickness  of  this  coating  was  estimated 
to  be  of  the  order  of  2  x  1G4  A  near  the  stagnation  point,  with  the  thick¬ 
ness  decreasing  with  distance  from  the  stagnation  point.  The  additional 
thermal  resistance  resulting  from  these  thin  coatings  was  negligible. 

The  data  for  any  single  probe  and  surface  were  not  all  obtained 
consecutively,  but  the  probe-surface  combinations  were  tested  in  a 
somewhat  random  pattern.  That  is,  each  probe  was  coated  and  cleaned 
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numerous  times  in  the  data  gathering  process.  The  data  repeatability 
was  quite  satisfactory  for  the  copper  surface.  For  the  noncatalytic 
surfaces,  there  was  some  scatter  in  the  data,  which  can  be  attributed 
to  the  inability  to  obtain  identical  surface  catalytic  properties  with 
successive  preparations  of  the  surfaces. 


SECTION  IV 

EXPERIMENTAL  RESULTS 


4.1  TREATMENT  OF  DATA 


The  most  common  method  of  presenting  experimental  heat  transfer 
data  at  low  Reynolds  numbers  is  to  normalize  the  results  with  the  first- 
order  boundary  layer  theory  applicable  at  high  Reynolds  numbers,  that 
is,  to  present  the  ratio  q/qBL,.  The  correlation  equation  of  Rosner 
(Ref.  47)  was  selected  for  the  calculation  of  q^L,: 


q  =  0.763 


(1) 


where  the  magnitude  of  the  enthalpy  difference  across  the  boundary 
layer  is 

Ah  =  (h0  -  h[tW)  —  (1  —  )  Ah0hem,  max  (2) 


This  equation  is  of  a  form  similar  to  the  correlation  equations  of  Fay 
and  Riddell  (Ref.  2)  and  Goulard  (Ref.  3}  for  the  frozen  boundary  layer, 
while  differing  from  Rosner’s  earlier  result  (Ref.  4)  only  by  the  inclusion 
of  the  factor  accounting  for  variable  pi u.  The  extent  of  recombination  <f> , 

defined  by  =  (3) 

has  been  generalized  by  Rosner  to  include  recombination  in  the  gas  phase 
as  well  as  at  the  body  surface.  For  a  completely  frozen  boundary  layer, 

$  can  be  reduced  to  the  simple  form 

*  =  —Hr-  (4) 

i  +  — 

k  w 

using  the  notation  of  Goulard  (Ref.  3).  The  constant  kw,  known  as  the 
catalytic  activity  or  surface  catalytic  reaction  rate  constant,  is  a  prop¬ 
erty  having  the  dimensions  of  velocity.  The  quantity  S,  which  is  a  meas¬ 
ure  of  the  diffusion  velocity  in  the  boundary  layer,  is  given  by 

S  =  0.763  (/9  pefie)°'5  (Sc)"”'*  pjl  (5) 


The  Prandtl,  Schmidt,  and  Lewis  numbers  were  taken  as  constants, 
since  the  uncertainty  of  their  correct  values  is  probably  greater  than 
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their  variation  across  the  boundary  layer.  The  Prandtl  number  was 
taken  as  0.71,  a  value  well  accepted  in  the  literature.  Following  Lees 
(Ref.  1),  the  Schmidt  number  was  taken  to  be  0.49.  The  Lewis  number, 
or  the  ratio  of  the  Prandtl  number  to  the  Schmidt  number,  was  then  cal¬ 
culated  to  be  1.45,  This  value  is  in  agreement  with  the  work  of  Brokaw 
(Ref.  48}.  (Many  writers  have  used  the  value  1. 4  for  the  Lewis  number, 
possibly  because  of  the  excellent  publicity  this  number  has  received 
from  other  sources.) 

The  inviscid  velocity  gradient  at  the  nose  was  calculated  using  the 
equation  of  Probstein  (Ref.  49)  for  modified  Newtonian  flow: 


Potter  and  Bailey  (Ref.  50)  have  found  that  the  modified  Newtonian  pres¬ 
sure  distribution  is  accurate,  even  at  these  low  Reynolds  numbers,  in 
the  region  O<0<  70  deg  for  highly  cooled  hemispheres. 

Since  Eq.  (1)  predicts  only  the  stagnation  point  heat  flux,  the  heat 
transfer  distribution  around  the  hemisphere  must  be  known  in  order  to 
predict  the  total  heat  transfer  to  the  hemisphere  for  comparison  with  the 
experimental  results.  The  heat  transfer  distribution  given  by  Lees 
(Ref.  1)  was  used  for  this  calculation.  The  assumptions  implicit  within 
Lees'  derivation  include  a  thin  boundary  layer  with  no  second-order 
effects,  cool  wall,  Prandtl  and  Lewis  numbers  of  unity,  constant  pju, 
and  chemical  equilibrium  at  every  point  in  the  boundary  layer.  Obviously, 
some  of  these  assumptions  are  inaccurate  for  the  present  experiments, 
and  the  resulting  effects  on  the  heat  transfer  distribution  must  be 
considered. 

Kemp,  Rose,  and  Detra  (Ref.  51),  using  the  "local  similarity"  con¬ 
cept  and  an  extension  of  the  ideas  used  by  Fay  and  Riddell  (Ref.  2), 
derived  another  expression  for  the  ratio  of  local  heat  transfer  rate  to 
stagnation  point  rate.  They  showed  that  the  effect  of  considering  pp 
variable  was  practically  equal  and  opposite  to  the  effect  of  using  0.  71 
instead  of  unity  for  the  Prandtl  number.  They  further  showed  that  the 
effect  of  non-unity  Lewis  number  was  small  for  Lewis  numbers  of  the 
order  of  1.4  and  that  the  effect  of  a  nonequilibrium  boundary  layer  was 
small  provided  the  wall  was  fully  catalytic.  Dorrance  (Ref.  52)  has  shown 
that  the  results  of  Kemp,  Rose,  and  Detra  are  equivalent  to  the  distri¬ 
bution  of  Lees  multiplied  by  a  factor  which  is  a  function  only  of  a  pres¬ 
sure  gradient  parameter.  It  is  further  shown  that  for  most  body  shapes 
of  practical  interest  (including  hemispheres)  this  factor  is  very  near 
unity.  The  shock  tube  data  of  Kemp,  Rose,  and  Detra  agree  equally  well 
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with  their  own  predictions  and  the  distribution  of  Lees.  In  addition,  Vidal 
and  Wittliff  (Ref.  53)  and  Hickman  (Ref.  54)  have  presented  data  at  low 
Reynolds  numbers  which  are  in  good  agreement  with  the  distribution  of 
Lees  for  hemispheres.  Their  results  show  that  the  heat  transfer  distri¬ 
bution  is  insensitive  to  Reynolds  number  for  hypersonic  flow. 

The  above  considerations  strongly  support  the  use  of  Lees’  heat 
transfer  distribution  in  the  treatment  of  the  present  data  for  a  fully 
catalytic  surface.  However,  for  the  " none ata lytic”  surfaces,  there  are 
no  such  firm  grounds  available  to  justify  its  use.  Nonetheless,  it  will  be 
used  for  these  surfaces  also,  since  no  better  information  exists. 

4.2  SURFACE  BEHAYIOR 

Before  the  experimental  results  are  presented,  it  is  appropriate 
to  comment  on  the  apparent  behavior  of  the  surfaces  when  exposed  to  the 
tunnel  flow.  The  change  of  surface  catalytic  properties  with  time  during 
the  test  runs  was  evident  to  some  degree  in  almost  all  cases  and  was 
of  particular  concern  with  the  coated  surfaces. 

The  test  procedure  generally  followed  was  to  begin  the  measure¬ 
ments  with  the  probe  located  at  the  station  farthest  downstream  from  the 
nozzle,  where  the  heating  rate  was  least.  After  this  measurement  was 
obtained,  the  probe  was  moved  from  station  to  station  toward  the  nozzle, 
until  the  final  station  at  the  nozzle  exit  plane  was  reached.  Then  this 
traverse  was  repeated,  usually  in  both  forward  and  reverse  directions. 

The  bare  copper  surface  exhibited  an  increase  of  about  15  percent 
in  total  heat  transfer  at  the  initial  station  after  one  complete  traverse  of 
the  test  section.  The  time  required  for  a  complete  traverse  of  the  test 
section  was  approximately  twenty  minutes.  After  this  first  traverse  of 
the  test  section,  there  was  no  observable  change  in  the  heat  transfer  rate 
at  a  particular  station  no  matter  how  many  additional  traverses  were 
made. 

The  fact  that  some  change  in  the  probe  surface  had  occurred  was 
evident  from  visual  inspections.  The  probe,  which  initially  had  a  shiny 
appearance  from  the  polishing  with  emery  paper,  had  distinctly  changed 
in  appearance  following  a  series  of  measurements  in  the  wind  tunnel. 

The  altered  surface  appeared  dark  and  dull  to  the  eye,  and  it  appeared 
that  this  new  surface  was  actually  a  thin  deposit  of  some  material.  A 
spectrographic  analysis  of  this  deposit  revealed  the  characteristic  lines 
of  copper,  possibly  in  an  oxide  form. 
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The  most  probable  source  of  this  copper  deposit  is  the  copper  anode 
in  the  arc  heater.  Because  of  the  action  of  the  arc  column,  the  anode 
does  erode  and  must  be  replaced  periodically,  say,  once  or  twice  a 
week.  This  erosion  is  so  slow  that  the  flow  contamination  is  completely 
~:gligible  for  most  aerodynamic  measurements.  On  the  other  hand,  when 
the  measurement  is  sensitive  to  minute  changes  in  the  properties  of  the 
body  surface,  as  in  catalytic  heat  transfer  measurements,  a  thin  deposit 
of  copper  could  significantly  alter  the  microscopic  surface  structure  of 
the  body  as  seen  by  the  gas  particles  themselves,  thereby  affecting  the 
catalytic  properties  at  the  wall. 

It  is  conjectured  that  in  the  present  experiments  a  very  thin  layer 
of  copper  was  deposited  on  the  probe  surface,  leaving  the  surface  rather 
porous  and  "rough11  on  the  microscopic  scale  as  compared  with  the  orig¬ 
inal  "smooth"  surface.  This  "microscopically  rough"  surface,  with  its 
peaks  and  crevices,  would  contain  a  greater  effective  surface  area  and 
would  have  the  ability  to  "adsorb"  gas  particles  for  longer  periods  of 
time,  thereby  enhancing  the  effective  catalytic  properties  of  the  surface. 
Since  copper  is  known  to  have  good  catalytic  surface  properties,  the 
properties  of  this  "roughened"  copper  surface  should  be  excellent  for 
catalytic  recombination  of  atoms. 

The  surfaces  coated  with  silicon  monoxide  and  ortho- phosphoric 
acid  exhibited  the  same  increase  in  heat  transfer  rate  with  time  as  did 
the  bare  copper  surface,  except  that  the  change  did  not  stop  after  the 
first  traverse  of  the  test  section,  but  continued  with  each  traverse  until 
the  heat  transfer  rates  to  the  bare  copper  were  approached.  This  degen¬ 
eration  of  the  noncatalytic  surface  was  apparently  caused  by  the  deposi¬ 
tion  of  copper  rather  than  the  erosion  of  the  surface  coating,  since  the 
electrical  insulating  property  of  the  silicon  monoxide  was  the  same  before 
and  after  each  run. 

Since  it  was  intended  for  the  bare  copper  probe  to  demonstrate  the 
effects  of  a  highly  catalytic  surface,  those  measurements  obtained  during 
the  initial  traverse  of  the  test  section  were  discarded.  Only  those  measure¬ 
ments  obtained  after  the  heat  transfer  rates  became  stable  are  included 
in  this  report.  Likewise,  the  purpose  of  the  coated  probes  was  to  provide 
measurements  for  a  highly  noncatalytic  surface.  Hence,  only  those  meas¬ 
urements  taken  during  the  first  traverse  of  the  test  section  are  included 
in  this  report.  Even  these  first  measurements  include  some  effect  of  the 
changing  surface  properties  caused  by  the  forming  layer  of  copper.  The 
extent  of  this  effect  will  be  discussed  when  the  data  are  presented. 

The  observed  variation  in  the  catalytic  properties  of  these  surfaces 
is  not  surprising  in  view  of  the  experience  reported  by  others.  In 
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particular,  Prok  (Refs.  55  and  56)  has  found  that  surface  catalytic  prop¬ 
erties  of  metals  and  metallic  oxides  are  strongly  dependent  on  the  type 
of  surface  preparation  as  well  as  the  duration  of  exposure  to  atomic 
specie.  He  has  also  found  that  the  thickness  of  an  evaporated  film  and 
the  type  of  substrate  were  quite  important  in  surface  catalysis.  Also, 
Myerson  (Ref.  57)  reported  that  the  catalytic  efficiency  of  silver  to  oxy¬ 
gen  atom  recombination  was  initially  a  direct  function  of  its  duration  of 
exposure  to  atomic  oxygen  fluxes.  Hartunian  and  Thompson  (Ref.  13) 
also  found  that  pre-exposure  of  their  probes  to  dissociated  flow  was 
necessary  to  produce  reliable  data. 

4.3  EFFECT  OF  FINITE  CATALYTICITY 

Although  a  number  of  independent  investigators  have  reported  values 
of  surface  catalytic  properties  in  the  literature,  Rosner  {Ref.  58)  has 
correctly  cautioned  against  the  use  of  these  recombination  coefficients 
for  the  composite  conditions  likely  to  be  met  in  aerodynamic  practice. 
Therefore,  the  recombination  coefficients  for  the  surfaces  used  in  the 
present  investigation  were  not  known  a  priori,  and,  indeed,  a  signifi¬ 
cant  portion  of  the  data  analysis  consisted  of  determining  the  proper 
values  to  be  used  in  comparing  the  experimental  results  with  theory. 

Thus  far,  two  parameters  which  can  be  related  to  wall  catalytic 
properties  have  been  introduced.  These  are  the  extent  of  recombination, 
ip,  and  the  surface  reaction  rate  constant,  kw,  related  by  Eq.  (4).  The 
limiting  case  of  "complete"  recombination,  <f>  =  1,  corresponds  to  an 
infinite  value  of  kw,  while  the  case  of  zero  wall  recombination,  =  0, 
corresponds  to  kw  =  0.  For  finite  (realistic)  values  of  kw,  the  heat 
transfer  rates  will  fall  between  these  two  limiting  cases.  In  order  to 
show  the  effect  of  finite  values  of  kw,  calculations  have  been  made 
using  Eqs.  (1)  through  (6)  with  the  present  flow  properties  for  kw  values 
of  102,  1()3,  10^,  and  105  cm/sec.  The  results  of  these  calculations  are 
compared  with  the  limiting  case  $  =  1  in  Fig.  5.  It  is  evident  that  a  sur¬ 
face  with  kw  near  102  cm/ sec  would  behave  essentially  as  a  completely 
noncatalytic  surface,  whereas  a  surface  with  kw  near  10^  cm/ sec  would 
act  essentially  as  a  "fully  catalytic"  surface.  Only  for  values  of  kw 
between  these  limits  would  the  effects  of  finite  atom  recombination  be 
detected. 

A  third  catalytic  parameter  which  is  of  interest  is  the  recombination 
coefficient,  yw,  defined  as  the  ratio  of  the  actual  frequency  of  surface 
encounters  leading  to  recombination  to  the  total  interfacial  collision  fre¬ 
quency  for  atoms.  In  other  words,  yw  is  the  probability  that  an  atom 
will  recombine  upon  striking  the  wall.  The  recombination  coefficient  is 
dimensionless  and  has  a  microscopic  interpretation,  as  opposed  to  the 
phenomenological  interpretation  of  the  surface  reaction  rate  constant,  kw . 
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With  the  assumption  of  a  Maxwellian  velocity  distribution  for  the 
atoms  striking  the  surface,  the  relation  between  yw  and  kw  has  been 
given  by  Goulard  (Ref.  3)  and  Rosner  (Ref.  59)  as 


Since  the  upper  limit  of  yw  is  unity,  there  will  be  a  finite  upper  limit  to 
kw ,  according  to  Eq.  (7),  for  each  wall  temperature.  Therefore,  the 
limiting  case  of  "complete  recombination",  with  0=1  and  kw  infinite, 
is  not  realistic  in  the  strict  sense.  This  can  be  seen  in  Pig.  5,  where  the 
heat  transfer  rates  for  several  values  of  yw  are  compared  with  the  case 
0=1.  The  effect  of  using  y w  =  1  instead  of  0  =  1  is  to  reduce  the  predicted 
heat  transfer  results  by  a  few  percent.  In  the  correlation  of  experimental 
measurements,  however,  the  difference  between  the  cases  0=1  and 
yw  =  1  may  not  be  of  great  significance  when  the  usual  precision  of  heat 
transfer  measurements  is  considered. 

4.4  RESULTS  FOR  COPPER  SURFACES 

The  measured  total  heat  transfer  rates  to  the  copper  surfaces  are 
compared  in  Fig.  6  with  the  values  calculated  from  Eq.  (1)  and  Lees' 
distribution.  The  calculations  were  made  assuming  "complete  recom¬ 
bination"  at  the  surface,  that  is,  0=1.  Also  presented  in  the  figure  are 
the  q/qgL  rati°s  predicted  by  the  second-order  theories  of  Cheng(Ref.  30) 
and  Van  Dyke  (Ref.  41).  The  location  of  the  curve  representing  Van  Dyke's 
theory  is  generally  typical  for  several  other  second-order  boundary  layer 
theories.  The  theory  of  Cheng  is  unique  in  that  it  predicts  a  higher  heat 
transfer  rate  than  most  other  theories;  yet  it  approaches  the  free-molecule 
limit  as  the  Reynolds  number  decreases. 

There  is  excellent  agreement  between  the  experimental  data  and  the 
theory  of  Cheng  under  the  simplifying  assumption  0  =  1.  The  indication  is 
that  the  bare  copper  surface,  when  exposed  to  the  tunnel  flow  long  enough 
for  the  transient  surface  effects  to  disappear,  is  an  extremely  efficient 
catalyst  to  nitrogen  atom  recombination. 

Since  it  was  stated  in  the  previous  section  that  the  more  realistic  case 
of  full  recombination  is  represented  by  yw  =  1  instead  of  0  =  1,  these 
same  experimental  data  are  compared  with  the  theory  for  the  case  yw  =  1 
in  Fig.  7.  Under  this  condition,  the  calculated  values  of  0  range  between 
0.  861,  at  the  lowest  Reynolds  number,  and  0.948  at  the  highest  Reynolds 
number.  The  effect  is  to  increase  the  values  of  q/^BL  slightly,  since 
qgL  would  be  decreased  according  to  Fig.  5. 
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The  data  presented  in  Fig.  7  for  the  case  7W  =  1  fall  slightly  on  the 
high  side  of  Cheng's  curve,  but  they  are  in  no  less  agreement  with  the 
theory  than  are  the  data  of  Fig.  6  for  $  =  1.  It  is  of  interest  to  note  that 
if  kw  were  calculated  from  the  expression  given  by  Cheng  (Ref.  30) 


instead  of  using  Eq.  (7),  then,  with  yw  =  1,  the  data  would  fall  approxi¬ 
mately  half-way  between  their  respective  locations  on  Figs.  6  and  7,  thus 
coinciding  with  Cheng's  curve. 

The  experimental  uncertainty  in  these  measurements  is  estimated 
to  be  about  ±5  percent,  mainly  because  of  the  probable  uncertainty  in  the 
total  enthalpy  measurement  and  the  use  of  Lees'  distribution.  The 
observable  experimental  scatter  is  also  of  this  same  order.  It  is  felt 
that  the  data  as  presented  in  Fig.  7  represent  the  best  comparison  with 
theory  since  this  figure  does  represent  the  most  realistic  catalytic  sur¬ 
face  condition.  Certainly  within  the  limits  of  experimental  accuracy,  the 
bare  copper  surface  behaved  as  a  fully  catalytic  surface,  with  7W  =  1. 

4.5  RESULTS  FOR  COATED  SURFACES 

The  measured  heat  transfer  rates  to  the  surfaces  coated  with  silicon 
monoxide  and  ortho -phosphoric  acid  are  presented  in  Fig.  8.  These 
results  have  been  normalized  with  the  values  predicted  by  multiplying 
the  results  of  Eq.  (1),  with  <j>  -  1,  by  the  values  of  q/qBL  obtained  from 
Cheng's  second-order  theory.  Also  presented  in  the  figure  are  lines 
representing  constant  values  of  kw  and  yw .  The  results  for  the  bare 
copper  surfaces,  if  presented  in  this  figure,  would  lie  between  the  top 
margin  of  the  figure  and  the  line  yw  =  1.0. 

There  is  a  marked  reduction  in  heat  transfer  with  the  coated  sur¬ 
faces,  being  as  much  as  30  to  40  percent  at  the  lowest  Reynolds  numbers. 
The  upward  trend  in  these  data  with  increasing  Reynolds  number,  evident 
in  Fig.  8,  would  be  expected  to  follow  the  slope  of  lines  representing  con¬ 
stant  values  of  yw  if  there  were  no  change  in  the  surface  catalytic  prop¬ 
erties  during  the  run.  The  amount  by  which  the  actual  slope  of  the  data 
from  a  particular  traverse  of  the  test  section,  exceeds  the  slope  of  a 
typical  line  representing  constant  yw,  say  yw  =  0. 1,  is  a  measure  of 
the  effect  of  the  changing  surface  properties  with  increasing  exposure 
time.  From  Fig.  8,  this  effect  is  responsible  for  about  ten  percent  of 
increase  in  heat  transfer,  during  one  test  section  traverse,  above  the 
increase  which  would  result  from  the  condition  of  constant  7W . 
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These  results  clearly  indicate  the  reduction  in  heat  transfer  which 
is  possible  in  the  frozen  boundary  layer  regime  when  a  noncatalytic  sur¬ 
face  is  employed.  Typically,  values  of  the  recombination  coefficient  -yw 
range  between  0.  01  and  0. 1,  compared  with  1.  0  for  the  catalytic  surfaces. 
The  large  scatter  in  the  data  results  from  the  inability  to  obtain  repeat- 
able  surface  properties  each  time  a  surface  was  prepared. 

It  is  interesting  to  note  that  the  initial  measurements  with  a  freshly 
polished  copper  probe,  taken  at  the  lowest  Reynolds  number  before  pro¬ 
longed  tunnel  exposure,  were  about  15  percent  below  the  final  steady- 
state  values,  corresponding  (from  Fig.  8)  to  a  value  of  yw  of  roughly  0.6. 
This  value  is  in  good  agreement  with  the  results  of  Ammann  (Ref.  60)  for 
nitrogen  atom  recombination  on  copper. 

4.6  COMPARISON  WITH  PREVIOUS  MEASUREMENTS 

Experimental  data  for  local  stagnation  point  heat  transfer  at  low 
Reynolds  numbers  have  been  reported  by  several  investigators  (Refs.  53, 
54  and  61  through  66).  As  previously  stated,  there  have  been  some  dis¬ 
crepancies  between  experiments  from  different  laboratories.  It  is 
hoped  that  the  present  results  will  aid  in  the  reconciliation  of  these 
discrepancies. 

Ferri  and  co-workers  (Ref.  62)  have  presented  experimental  data 
for  the  vorticity  interaction  regime  for  stagnation  temperatures  of 
1600  and  2300  °R.  These  data,  when  presented  in  the  form 
show  an  increasing  vorticity  effect  for  increasing  stagnation  tempera¬ 
ture.  These  data  are  in  general  agreement  with  the  theory  of  Ref.  62 
and  also  with  the  theory  of  Cheng,  with  both  experiment  and  theory 
displaying  a  vorticity  effect  at  much  higher  Reynolds  numbers  than 
predicted  by  other  second-order  theories. 

Other  measurements  for  the  vorticity  interaction  regime  have  been 
reported  by  Hickman  (Ref.  54).  These  measurements,  obtained  at  a 
stagnation  temperature  corresponding  to  room  temperature,  are  in 
reasonable  agreement  with  the  theory  of  Van  Dyke  and,  therefore,  indi¬ 
cate  a  much  smaller  vorticity  effect  than  was  observed  by  Ferri  et  al. 
Even  though  the  data  of  Hickman  and  Ferri  et  al.  do  not  agree,  they  are 
consistent  in  the  sense  that  they  show  an  increasing  vorticity  effect  with 
increasing  stagnation  temperature. 

As  the  Reynolds  number  is  further  decreased  below  those  corre¬ 
sponding  to  the  vorticity  interaction  regime,  viscous  effects  in  the  layer 
between  the  shock  and  the  body  become  important.  Ferri  and  Zakkay 
(Ref.  63)  have  also  reported  experimental  measurements  in  this  viscous 
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layer  regime  for  a  stagnation  temperature  of  2100°R.  These  results  are 
in  agreement  with  the  earlier  data  of  Ferri  et  al.  and  show  a  decrease  in 
q /  qj3L  in  accord  with  the  prediction  of  Cheng. 

Potter  and  Miller  (Ref.  66)  have  also  reported  measurements 
obtained  earlier  in  the  same  wind  tunnel  used  in  the  present  measure¬ 
ments.  Their  measurements  were  also  obtained  in  nitrogen,  but  at  a 
higher  stagnation  pressure  so  that  the  flow  contained  no  dissociation. 
These  data  are  in  good  agreement  with  the  Cheng  theory  presented  in 
the  normalized  form  as  in  Fig.  9. 

Shock  tunnel  measurements  at  low  Reynolds  numbers  have  been 
reported  by  Vidal  and  Wittliff  (Ref.  53).  There  is  considerable  scatter 
in  these  measurements  and  they  could  be  used  equally  to  support  all 
theoretical  predictions  in  the  vorticity  interaction  regime.  At  the  lowest 
Reynolds  numbers,  the  data  are  in  general  agreement  with  Cheng's 
theory  when  the  enthalpy  associated  with  the  free-stream  atom  concen¬ 
tration  is  neglected.  Although  the  models  used  were  presumed  to  have 
completely  noncatalytic  surfaces,  some  of  the  observed  scatter  could 
be  caused  by  a  variation  in  surface  catalytic  properties  similar  to  that 
observed  in  the  present  investigation. 

The  data  of  Ferri  et  al.  (Refs.  62  and  63)  and  Potter  and  Miller 
(Ref.  66)  are  reproduced  in  Fig.  9.  These  data  may  be  compared  with 
the  present  results  given  in  Fig.  7.  Agreement  between  the  data  in  these 
two  figures  is  relatively  good. 


SECTION  V 
CONCLUSIONS 


This  report  has  presented  the  details  and  results  of  an  experimental 
investigation  of  heat  transfer  in  nonequilibrium  dissociated  hypersonic 
flow  with  surface  catalysis  and  second- order  boundary  layer  effects. 

This  section  is  devoted  to  an  enumeration  of  the  significant  observations 
made  during  the  investigation  and  conclusions  drawn  from  the  investigation. 

Since  a  significant  effort  was  devoted  to  calibration  of  the  nozzle 
flow,  two  observations  concerning  the  flow  diagnosis  will  be  made.  First, 
the  combination  mass-flow  and  total- enthalpy  probe,  although  not  entirely 
free  from  error,  provided  a  satisfactory  means  of  determining  locally 
the  total  enthalpy  of  the  flow.  Secondly,  the  nitric  oxide  titration  tech¬ 
nique  proved  to  be  a  useful  method  of  determining  the  free-stream  atom 
concentration,  even  in  high  speed  flows. 
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The  major  observations  and  conclusions  from  the  results  of  the 
experimental  heat  transfer  measurements  are  summarized  as  follows: 

1.  The  bare  copper  surfaces  exhibit  an  initial  recombination 
coefficient  of  about  0.6.  The  surface  recombination  increases 
rapidly  with  time  until  no  further  time-dependence  is  observed. 
At  this  steady-state  condition,  surface  recombination  is  virtu¬ 
ally  complete,  that  is,  yw  =  1.  This  transient  effect  is  attributed 
to  an  enhancing  of  the  microscopic  surface  catalytic  properties 
owing  to  the  formation  of  a  fine  deposit  of  copper  on  the  probe 
surface. 

2.  The  experimental  results  for  copper  surfaces,  normalized  by 
the  values  computed  from  Rosner's  stagnation  point  heat  transfer 
equation  assuming  yw  =  1,  are  in  agreement  with  the  thin  shock 
theory  of  Cheng.  Also,  they  are  in  agreement  with  the  experi¬ 
mental  measurements  of  Perri  and  Zakkay. 

3.  A  significant  reduction  in  heat  transfer  to  a  noncatalytic  surface 
has  been  demonstrated.  The  coated  surfaces  used  in  this  inves¬ 
tigation  were  highly  susceptible  to  catalytic  property  changes 
with  time  of  exposure. 

4.  The  use  of  differential  catalytic  probes  for  determining  free- 
stream  atom  concentrations  is  currently  under  study  at  several 
laboratories.  The  transient  surface  effects  observed  in  this 
investigation  point  out  a  potentially  serious  limitation  in  the  use 
of  such  probes  to  measure  absolute  values  of  atom  concentra¬ 
tion.  To  overcome  this  difficulty,  surfaces  must  be  used  which 
retain  known  catalytic  properties  over  long  exposure  periods. 

5.  The  results  of  this  investigation  have  by  no  means  ended  the 
need  for  additional  experimental  study  of  heat  transfer  with 
surface  atom  recombination  and  low  density  effects.  To  the 
writer's  knowledge,  this  effort  represents  the  first  instance 
where  combined  surface  atom  recombination  and  second-order 
boundary  layer  effects  have  been  studied  experimentally. 
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Fig.  2  Toto I- Enthalpy  and  Mass-Flow  Probe  Installed  in  Tunnel 
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Fig*  4  Water-Cooled  Heat  Transfer  Probe 
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Fig.  6  Heat  Transfer  to  Copper  Hemispheres  for  <ft 
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